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)Q ■ ABSTRACT 

\0 • 

This work is part of the NUGA survey of CO emission in nearby active galaxies. We present observations of NGC 4569, a member of the 
lO Virgo Cluster We analyse the molecular gas distribution and kinematics in the central region and we investigate a possible link to the strong 
starburst present at the nucleus. 70% of the [1.1 ± 0.2] x 10' Mq of molecular gas detected in the inner 20" is found to be concentrated within 
1 the inner 800 pc and is distributed along the large scale stellar bar seen in near-infrared observations. A hole in the CO distribution coincides 
' with the nucleus where most of the Ha emission and blue light are emitted. The kinematics are modelled in three different ways, ranging 
^ , from the purely geometrical to the most physical. This approach allows us to constrain progressively the physical properties of the galaxy and 
T ""j eventually to emerge with a reasonable fit to an analytical model of orbits in a barred potential. Fitting an axisymmetric model shows that the 
I^nJ non-circular motions must be comparable in amplitude to the circular motions (120km s"'). Fitting a model based on elliptical orbits allows us 
to identify with confidence the single inner Lindblad resonance (ILR) of the large scale bar. Finally, a model based on analytical solutions for 
the gas particle orbits in a weakly barred potential constrained by the ILR radius reproduces the observations well. The mass inflow rate is then 
estimated and discussed based on the best fit model solution. The gravitational torques implied by this model are able to efficiently funnel the 
gas inside the ILR down to 300 pc, although another mechanism must take over to fuel the nuclear starburst inside 100 pc. 

Key words, galaxies: individual: NGC 4569 - galaxies: active ~ Galaxy: kinematics and dynamics 



1. Introduction 



NUGA, for NUclei of GAIaxies dGarcfa-Burillo et aLll2003al) . 
is a survey of nearby active galaxies aiming for a deep insight 
into the mechanisms responsible for feeding the nuclei of ac- 
tive galaxies. High resolution (down to 0.5") and high sensitiv- 
ity (rms of ~2mJy beams"' for lOkms"' channels) observa- 
tions of 12 galaxies were obtained with the IRAM Plateau de 
Bure Interferometer (PdBI) and the IRAM 30 m telescope. The 
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* Based on observations carried out with the IRAM Plateau de Bure 
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(Germany) and IGN (Spain). 



CO 7=1-0 and 2-1 transitions thus mapped allow us to study 
the dynamics of the molecular gas in the central kiloparsec of 
the galaxies with a linear resolution as high as ~10pc. 

The main results achieved by the NUGA project thus far 
can be summarized as follows: 

- A great variety of mor phologies and dynamics is ob served: 



2 and 1-arm spirals dGarcfa-Burillo et al J [2003b), bars 



rings ( Combes et al. 20041). circum nuclear discs, asymme 



tries, and warps dKrips et alj|200l . As a corollary, there is 
no evidence for features (at scales >10pc) that would be 
typical of active galaxies and that could be uniquely linked 
to the activity of the nucleus. 
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- We do not see any 2D-kinematic pattern indicative of 
systematic inflow motions in the circumnuclear disks of 
some of the NUGA galaxies examined thus far down 
to the scales of our spatial resolution (10 pc at best). In 
some cases, e.g NGC4826 (iGarcia-Burillo et all l2003bl) 



and NGC7217 (ICombes et al 



20041) . the observed pertur- 



bations to the disks seem less hkely to drive than to inhibit 
the fueling of the nucleus. 
- A detailed study of the gravitational torques exerted by the 
large scale disk asymmetries on the gas in four galaxies 
showed that such asymmetries may not be the dominant 
mechanism for angular momentum removal in the nuclear 
regions (jGarcfa-Burillo et al. 2005). 

These first results suggest either that low luminosity AGN 
(the NUGA sample includes Seyferts and LINERs) do not re- 
quire efficient fueling to sustain their luminosity (see e.g. IHoI 
20031) . or that the timescales for fueling and the onset of activ- 
ity are so different that b oth cannot be observed simultaneously 
(see e.g. ICombesI 120041) . However, more galaxies in the sam- 
ple need to be analysed before any definitive conclusions about 
AGN fueling can be drawn. 

NGC4569 (M 90) is a bright SAB(r s )ab galaxy a t a dis- 



1983 



tance of 17Mpc jBottineUi et al 
Virgo Cluster. A large s cale bar is seen in 



Tullvl 119881) in the 
NIR 



images 

(iLaurikainen & Saloll2002h and is almost aligned with the ma- 



jor ax is of the galaxy (R^=15deg according to iJogee et al 



20051) . The galaxy harbour s a nucleus of the transition type 



(type T2 in iHo et al.l Il997h which exhibits, by far, the most 



pronounced nuclear starburst activity amo ng the LINERs an d 



transition nuclei with available UV data (IMaoz et al.llI998h . 
As such it became a prime target for studies of the ori- 
gin of the ionization in transition nuclei. The p resence of a 
supergiant-dominated s tarburst was established bv lKeel (Il996l) . 
Barth & Shieldsl dlOOO") demonstrated that a starburst with a 
large population of very hot Wolf-Rayet stars could produce 
the observed spec trum of the nucleus and AGN activity is 
ruled out by X- ray ( Terashima et alj |2000: Tschoke et al. 2001; 
Ho et al.1 200 lb and rad i o observations (Hummel et al.. ,1987 ; 
Neff & Hutchingslll992l) . iGabel & Bruhweilej ( |2002|) showed 
that for the starburst to be the only source of photoionization 
(no AGN) extreme conditions are required, with ~5xlO'* O and 

B stars squeezed int o the inner 30 p c. 

On larger scales iTschoke et al.l ( 1200 lb detected a large soft 



X-ray emitting region above the disk of NGC 4569 wh i ch im- 
plies a large energy input into the halo. iHensler et alJ (l2003h 
reported on a bipolar outflow seen in X-rays and a giant out- 
flow (~ 10 kpc in length) to the western side of the galaxy in Ho- 
(see also Bomans et al 2007, submitted). Using the Effelsberi 
radio telescope at 4.85 GHz and 8.35 GHz. lChyzy et alJ (l200 
discovered large symmetric lobes of polarized radio emission 
extending up to 24 kpc from the galactic disk on each side of 
the nucleus. These lobes may have been powered by the nuclear 
starburst. 

H I line observations reveal that the galaxy is H i deficient 
compared to galaxies of simila r type and it must have lost 
more than 90% of its atomic gas ( Warmelslll988 : Cavatte et al 
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Fig. 4. First moment map of the CO(l-O) emission. Regions 
with low S/N are masked out. The contour levels are spaced 
by 10km s ' . The color scale is given on the right hand-side in 
kms"'. 



1990l) . This was most likely stripped due to rapid motion of the 



galaxy through the intrac luster medium d VoUmer et al.l 12004 : 
isoseni et al]|2006l) . 

Interferom etric CO o b servat io ns of NGC 4569 were 
presented in iHelfer et al.l (120031) . I Jo gee et al.l (l2005h and 
akanishi et al. I (l2005h . In this article, higher resolution and 
higher sensitivity CO(l-O) and, for the first time, CO(2-l) line 
observations are presented, and the analysis is focused on in- 
vestigating a possible link between the molecular gas kinemat- 
ics in the inner kiloparsec and the nuclear starburst. The ob- 
servations are presented in Section|2] The gas kinematics are 
modeled based on analytical expressions for the gas orbits. 
Three dififerent models are compared to the observations start- 
ing with a pure geometrical model and proceeding to more re- 
alistic models in Section[3] Comparing the different models to 
the observations allows us to progressively constrain the kine- 
matic properties of the galaxy and emerge with a reasonable fit. 
The mass inflow due to gravitational torques based on the best 
fit version of the most realistic model is estimated in Section]?] 
The goal of this estimate is to quantify the possible contribution 
of the gravitational torques to the fueling of the nuclear star- 
burst. The results are discussed in Section]5]and summarized in 
Section]6] 

2. CO Observations and global properties of the 
emission 

2.1. Observations 

The galaxy was observed with the IRAM Plateau de Bure 
Interferometer (PdBI) in 2003 (B and C configurations) and 
2004 (A configuration) and with the IRAM 30 m telescope in 
2004 in the CO(l-O) (115 GHz) and CO(2-l) (230GHz) emis- 
sion lines. The PdBI and 30 m telescope receiver characteris- 
tics, the observing setup and the observing procedures are the 
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Fig. 1. a) CO(l-O) (grey scale) and CO(2-l) (contours) integrated maps with a clip at zero and no correction for primary beam 
attenuation. The color scale is given on the right hand-side in Jy beam"' km s"' . The contours start at 4 cr and are spaced by 10 cr 
with cr=0.38 Jy beam ' km s ' . The 115 GHz and 230 GHz beams are shown on the bottom right and left respectively. The axes 
give the angular offset from the phase center (symbolized by the cross) in arcseconds. We recall that, at 17 Mpc, 1" corresponds 
to 82 pc along the major axis, b) CO(2-1)/CO(1-0) ratio map with the CO(l-O) contours overlaid. The contours start at 4 cr and 
are spaced by 10 cr, where cr=0.27 Jy beam ' kms c) HST Ha greyscale with CO(l-O) contours. 



same as in lCombes et al.l (12004)) . The quasar 3C273 was used 
for PdBI gain calibration. The maps obtained from the 30 m ob- 
servations are used to compute the short spacings and complete 
the interferometric measurements u sing the SHO RT-SPACE 
task in the GILDAS software (Guillo teau & Lucasl BoOO) as in 
Combes et aD (|2004 . The RMS in 10km s ' wide channels is 
3.9 and 5.4 mJy beam ' for natural weighting and the beam- 
sizes are 2.5"xl.6" and l.l"x0.7" at 115 GHz and 230GHz, 
respectively. 



2.2. Morphology and mass of molecular gas 

Both CO(l-O) and CO(2-l) emission are distributed over a 
region 17"x6" in size (Fig. 1 a&b) roughly aligned with the 
major axis of the galaxy (PA=24deg). The CO(l-O) distribu- 
tion shows three peaks, one peak close to the center and two 
symmetrical peaks at ~5" from the center and aligned with 
the overall elongation. The south-western peak is the strongest. 
The central peak does not coincide with the dynamical center: it 
appears to be the top of a slightly curved ridge encircling a hole 
in the CO distribution centered ~1" northwest of the dynami- 
cal center Since our ob servations have a S/N ^ 10 times higher 
than those presented by lNakanishi et alJ (120051) . a weaker con- 
tinuous distribution is detected in addition to the three main 
peaks, even when going to the highest resolution in the CO(2- 
1) map. 

The integrated CO(l-O) flux within a radius of 20" is 
450Jykms ' with an uncertainty due to calibration errors 
of 15%. This is com patibl e with the fl uxes measured by 

(l2005h . The corre- 
t 0.2] X 10^ Mq using a conversion 
"2[Kkms ']"' dSolomon & Barrett 



Nakanishi et al.l (120051) and iJogee et al 



sponding H2 mass is [1.1 
factor ofX = 2.2 xl02°cm 



11991 ). 70% of this flux is concentrated within a radius of 10" 
(-800 pc). 



2.3. Velocity distribution 

The line-of-sight CO velocities span a range of 400 kms '. 
Channel maps are shown in Figures|2] and |3] The most no- 
ticeable features are: (1) the distribution is fairly symmetric 
with respect to systemic velocity and dynamical center; (2) the 
emission is elongated in most of the channels (from -150 to 
+ 150 kms"') along a direction close to that of the emission in 
the integrated map; (3) from one channel to the next as velocity 
increases, the elongated emission is slightly shifted toward the 
northwest; (4) the highest velocities (±200 kms"') are reached 
at ~5" from the nucleus, i.e. close to the two outer peaks in the 
integrated map; and (5) some low-level emission, mostly seen 
in the CO(l-O) channel maps (Fig.|2|l outside the elongated 
core, seems to follow the "spider" pattern of rotating disks, i.e. 
two "arms" on each side of the major axis shift from the north 
at low velocities to the south at high velocities. To guide the 
eyes of the reader the two arms are emphasized with thick red 
lines drawn by hand and overlaid to the channel maps in the 
bottom panel of Fig.|2] Due to the low level of the signal these 
arms appear discontinuous and the spider pattern is not seen 
in the velocity field of Fig. |4] (the first moment mixes this low 
signal with the noise). 

The fact that in individual channels the emission is elon- 
gated in a direction close to that of the major axis (PA=24deg) 
implies that significant amounts of gas must lie at "forbidden" 
velocities, i.e. velocities of a gi ven sign appear on both sides of 
the minor axis, as was noted by Jogee et al. ( 20051) . The veloc- 
ity field is shown in Fig. 13 it is clearly distorted compared to a 
disk with only circular rotation. In the central region a steep ve- 
locity gradient is running in a SE to NW direction, i.e. orthog- 
onal to the major axis. The kinematic line of nodes is therefore 
orthogonal to the major axis. This steep gradient results from 
the continuous shift of the emission in the channel maps (fea- 
ture (3) described above). 
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Fig. 2. CO(l-O) channel maps. The velocity relative to the systemic velocity of the galaxy (vhei=235kms" ) is given at the top left 
of each map. The phase center symbolized by the cross is at a'72000=12''36'"49.8*, 5y200()=13 °09'46.3". The beam is 1.8"xl.l" 
and the rms is 2.1 mJybeam The contours start at 4 cr and the levels are separated by 20 cr. In the bottom panel the channel 
maps are overlaid with thick lines emphasizing the two "arms" seen at low level (see Section l273T l. 
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Fig.3. CO(2-l) channel maps. The beam is l.l"x0.7" and the rms is 3.7 mJy beam '. The contours start at 4cr and the levels 
are separated by 20 cr. 

2.4. Line ratio a kinematic shock resulting from the asymmetrically perturbed 

gravitational potential. 



The CO(2-1)/(1-0) line ratio was computed after tapering the 
230 GHz data and restoring it with the same beam as for the 
115 GHz data; the ratio map is shown Fig.[T]b. The values are 
in the range 0.3-1.1 in temperature units. The highest values 
are reached in the ridge between the central and the southern 
peaks near the hole in CO. A high line ratio corresponds to 
high density and/or high temperature gas. A possible interpre- 
tation is that the molecular gas in the ridge is being compressed 
and heated by the nuclear starburst. This interpretation would 
also be consistent with the presence of a hole in the molecu- 
lar gas distribution: part of the molecular gas undergoing this 
pressure might have been dissociated or pushed out of the disk. 
A fraction of the initial molecular gas in the center must also 
have been consumed by star formation. The high ratio in the 
southern peak could be either due to a star forming region or 



2.5. Comparison to optical observations 

A composite color map from B, V and H broad band im- 
ages and the V-I color map from HST observations are pre- 
sented in Fig.|5] Dust lanes are clearly seen along the bar, 
which in NIR observation s extends up to a radius of ~5kpc 
(iLaurikainen & Saioll2002h . The bar appears to be aligned with 
the major axis of t he galaxy, like the CO emission, as noted by 
Jogee et al.l (120051) . Also striking is the difference in extinction 
between the eastern and western sides of the galaxy. This con- 
firms that the near side, i.e. the side that appears as the most 
extinguished because of a low contribution from foreground 
bulge stars, is the western side. This geometry is consistent 
with the spiral arms being trailing as in most spiral galaxies. 
The two outer CO(l-O) peaks coincide with peaks in dust ex- 
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Fig. 5. Left: RGB color image synthesized from B,V and H broad ban d images obtained from Observatoire de Haute Provence 
(OHP) and retrieved from the Goldmine database ( Gavazzi et alj|200"3 ). North is up and East is left. Right: V-I color index from 
HST WFPC2-camera observations (F555W and F814W filters, proposal 05375). The color scale is such that high V-I index 
appears in yellow and low V-I in blue. 



tinction, the southern one being the more prominent both in 
extinction and CO(l-O) emission. The central CO(l-O) peak 
does not coincide with any extinction peak, and the hole in 
the CO(l-O) distribution close to the location of the nucleus 
has blue optical colors. This agreement suggests that most of 
the starburst is concentrated within a region of ~100pc in ra- 
dius, supporting the hypothesis that the hole in the CO(l-O) 
emission is a consequence of this concentrated starburst hav- 
ing consumed or expelled the molecular gas from the center. 
The latter interpretation is further supported by an HST Ha im- 
age jPogge et al. I I2OOOI) which shows strong emission south of 
the hole and several features elongated toward the east (Fig.[T]i. 
Some Ha emission is expected due to the ionizing stars of the 
starburst but part of it might also trace the impact of the star- 
burst on the ISM and the elongated features to the east are sug- 
gestive of an outflow perpendicular to the plane. The east side 
being the far side, the counter-outflow (i.e. emerging from the 
other side of the plane) might be obscured and therefore not 
visible in Ha, although again visible further to the west as re- 
ported bv lHensler et al.l (l2003h . 



3. Kinematic modeling 

3.1. Motivation and strategy 

How can the CO emission described in the previous section 
be interpreted in terms of the radial distribution and kinemat- 



ics of the molecular gas in the plane of the galaxy? The two 
outer peaks seen in the integrated map could be due to a ring 
or to concentrations at the ends of a nuclear bar. The elon- 
gated shape seen in the channel maps could result from radial 
motions (streaming motions inside a density wave or due to 
shocks at the ends of a bar) or to elongated orbits. To discrimi- 
nate among these possible interpretations and obtain a quan- 
titative understanding of the deprojected gaseous disk, three 
different kinematic models are compared to the observations. 
The approach followed is progressive, starting with a pure geo- 
metrical model and proceeding to more 'physical' models, i.e. 
with more physical underlying hypotheses used to compute the 
molecular gas motions. This approach allows us to constrain 
the physical properties of the galaxy step by step and thus to 
converge more efficiently to a reasonable solution. This is more 
efficient because the most physical model is very sensitive to 
variations of the parameters and therefore more difficult to fit 
without a priori constraints. 



3.2. General description of tlie modeling 

The first model is based on circular orbits with additional radial 
motions, while the second model is constructed from elliptical 
orbits. These are purely geometrical hypotheses, and reality is 
expected to be more complex. However, the symmetry of the 
CO emission in the datacube (i.e. with respect to the dynamical 




Fig. 6. Each row shows a 50"x50" channel map of the observed data, the axisymmetric model, the ellipse orbit model and the 
barred potential model respectively (from left to right) at a given velocity. From top to bottom the velocities are -80, -50, 0, 50 
and 80kms"^. We recall that, at 17Mpc, 1" corresponds to 82 pc along the major axis. 




Fig. 7. Integrated maps (first row) and first moment maps (second row) for a) the observed data, b) the axisymmetric model, c) 
the elHpse orbit model and d) the barred potential model. In the observations, the velocity field can only be computed where the 
signal-to-noise ratio is high enough, i.e. within a chevron shaped field. For comparison we used the same mask when computing 
the velocity fields from the models. 
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Fig. 8. Deprojected integrated CO(l-O) emission of a) the axisymmetric model, b) the ellipse orbit model and c) the barred 
potential model. We recall that, at 17 Mpc, l"~82pc along the major axis. 



center and the systemic velocity) suggests the gas flow is well 
organized and dominated by ot=2 modes. Thus, if we assume 
a pure m=2 symmetry, comparing the observations to the first 
model can help to quantify the amplitude of non-circular mo- 
tions and estimate to what extent they can be reproduced by 
pure radial motions. Comparison to the second model allows 
us to evaluate the alternative interpretation of non-circular mo- 



tions, namely the possibility that the clouds populate near-to- 
closed orbits that are elongated. The third model is based on 
analytical solutions for gas particle orbits in a barred potential. 

The first two models are entirely specified by a radially 
varying set of pro perties similar to th at of the widely used 
'tilted ring model' ( Rogstad et alJll9'74 ). The values of the pa- 
rameters are specified at some galactocentric radii and linearly 
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Fig. 9. Orbit patterns in the inner 1.5 kpc of the elUptical orbit model (left) and the barred potential model (middle). Right-hand 
side: Resonance diagram, with the rotation frequency (Q) plotted in red, Q - k/2 in blue, and Q + k/2 in cyan. The rotation curve 
is overplotted in black with the scale on the right-hand axis. 



interpolated between those radii. The number of radii where 
the parameter values are specified and their values can be cho- 
sen freely, giving us the ability to reproduce any radial distri- 
bution to a good approximation. Simulated data are produced 
by randomly assigning gas particles to the disk according to 
the distributions specified, projecting them in the datacube and 
convolving by the synthesized beam. As most of the emission 
comes from within the PdBI primary beam full width at half 
maximum (43" at 1 15 GHz), the effect of primary beam atten- 
uation is neglected. For more details about the fitting procedure 
see Appendix lA. 1 1 

3.3. Axisymmetric model 
3.3.1. Results 

The parameters of the axisymmetric model (see details in 
Appendix lA. 21 ) reproducing the main features of the observa- 
tions with 7 radii are listed in Table[T] The corresponding chan- 
nel maps for 5 velocities are shown in the second column of 
Fig.|6] while the corresponding integrated intensity and veloc- 
ity maps are shown in Fig.|7]b. A face-on view of the model is 
shown in Fig.[8]a. 

The three peaks in the integrated map are reproduced by a 
ring with mean galactocentric radius of ~600pc and a central 
concentration. The column density must be higher by a factor 
of ~10 inside the ring (at 600pc) with respect to the density 
at 300 pc and 700 pc. The column density must increase again 
inside a radius of 300 pc to reach a value higher by a factor of 
~15 at the center The distribution of the CO emission in the 
channel maps requires an inward radial velocity of 120 km s"' 
in the ring at 600 pc. Lower radial velocities in the ring fail to 
reproduce the data. 

This model can reproduce the shape of the CO emission 
in the channel maps at all velocities. The slight shift of the 
emission toward the northwest from channel to channel and, 
as a consequence, the twisted isovelocity lines are also well 
reproduced (compare Fig.|7^ withlTJ)). Moreover, the ring in- 



terpretation provides a natural explanation for the fact that the 
elongated shape is aligned with the major axis of the galaxy: it 
is merely an inclination effect. The width of the emission along 
the minor axis, however, is larger than in the observations, and 
the two symmetrical peaks are not as highly contrasted as ob- 
served. Also, the radial velocity in the ring required to repro- 
duce the observations is extremely high and cannot be inter- 
preted as the manifestation of real radial motions. 

3.3.2. Interpretation 

We conclude from this first modeling that, although a ring of 
600 pc radius with a central concentration can explain to some 
extent the spatial distribution of the CO emission, the amplitude 
of the non-circular motions is unrealistically high, i.e. compa- 
rable to the amplitude of the rotation motions. No physics could 
explain the collapse of a ring at this speed . This result implies 
that pure radial motions alone cannot explain the departures 
from circular rotation, which must be at least partly due to the 
shapes of the orbits themselves. Paradoxically, while the best 
fit solution of the circular model reproduces the observations 
well, it implies that the model is not a good representation of 
reality and that elongated orbits must be considered. 



3.4. Elliptical orbit model 
3.4.1. Description 

Although elongated orbits are not expected to be elliptical in 
general, ellipses provide a first order approximation to any non- 
circular closed or nearly-closed shape. To estimate a possible 
contribution from the orbit elongation to the non-circular mo- 
tions and to get a picture of the orientations of the orbits in the 
galactic plane (the orbit pattern) we fit a model in which the gas 
particles populate elliptical orbits and their angular momentum 
is constant over the orbit. "Nearly-closed" orbits means orbits 
that are stable on the time scale of a rotation period, which 
is possible only if the orbits do not intersect, since otherwise 
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Table 1. Circular model parameters 



radius (pc) 


1.0 


170 


340 


560 


800 


1700 


2800 


rotation vel. (km s"') 


40.0 


30.0 


40.0 


50. 


70.0 


100.0 


110.0 


radial vel. (kms"') 


0.0 


0.0 


0.0 


-120.0 


0.0 


0.0 


0.0 


col. dens, (no units) 


6.0 


1.0 


0.5 


4.0 


0.4 


0.4 


0.4 


velocity disp. (kms"') 


40.0 


10.0 


10.0 


20.0 


10.0 


10.0 


10.0 


scale height (pc) 


20.0 


15.0 


10.0 


10.0 


10.0 


10.0 


10.0 



Table 2. Elliptical orbit model parameters 



charact. radius (pc) 


1.0 


85 


170 


510 


620 


680 


850 


1700 


2800 


5500 


charact. rot. vel. (kms"') 


40 


25 


10 


20 


40 


90 


120 


150 


200 


220 


ellipticity 


0.8 


0.8 


0.8 


0.7 


0.5 


0.3 


0.3 


0.3 


0.2 


0.2 


position angle (deg) 


110 


130 


160 


210 


160 


100 


20 


-180 


-230 


-330 


col. dens, (no units) 


3.0 


2.0 


0.2 


0.7 


0.8 


0.3 


0.03 


0.03 


0.03 


0.02 


velocity disp. (km s"') 


30 


35 


40 


40 


40 


30 


10 


10 


10 


10 


scale height (pc) 


20 


18 


17 


15 


10 


10 


10 


10 


10 


10 



the cloud shocks would re-orient the gas on different orbits. 
Thus, for this model to be self-consistent, the variations in po- 
sition angle and ellipticity must satisfy this 'no crossing' con- 
dition. The density of clouds on a given orbit is not uniform. 
If it were assumed independent of the rotation sense it would 
be inversely proportional to the velocity. However, it is known 
that molecular clouds are formed in density waves or shocks, 
from which stars form that then destroy the parent cloud. The 
density is therefore expected to increase abruptly when enter- 
ing the density wave or shock and to decrease smoothly after 
In this simple model, density waves correspond to overdensi- 
ties due to orbit crowding, i.e. are kinematic waves: they form 
at the aphelia of the ellipses when the position angles change 
with radius. To mimic the formation/destruction of the clouds 
we apply an exponential taper to the density. The application 
of this taper significantly improves the model fit. More details 
about this model are given in Appendix lA.3l 

3.4.2. Results 

The parameters of such a model best reproducing the obser- 
vations with 10 radii are listed in Table|2] The corresponding 
channel maps for 5 velocities are shown in the third column 
of Fig.|6] while the corresponding integrated map and veloc- 
ity field are shown in Fig.|2]c. The face-on view of the model 
solution is shown in Fig.[8]b and the orbital pattern in Fig.|9]a. 
Strong variations in the column density, the angular momen- 
tum, the orbit orientations and the velocity dispersion are re- 
quired at a radius of ~600pc to reproduce the peaks in the 
integrated map and the shape of the distribution in the chan- 
nel maps. As in the previous model the column density must 
increase by an order of magnitude at ~600pc, however inside 
the ~600pc radius the column density only drops by a factor 
of 4. A central concentration is still required within the inner- 
most 150pc. The angular momentum must decrease abruptly 
inside ~600pc. The velocity dispersion must be ~40kms ' at 
~600pc. The most remarkable result is that to reproduce the 
elongated core in the channel maps, it is necessary to invert 
the variation of the orbit position angles at the same radius of 



the column density peak (in the range 500-600 pc). The orbit 
position angle must be maximum at this radius and equal to 
~30 mod 180 deg. This feature of the orbit pattern seems to be 
strongly constrained by the observations. 

From Fig.|6] it can be seen that the low level emission is 
well reproduced by this model, giving us confidence that this 
emission is real and that the rotation velocity of the disk at 
5.5 kpc is ~22Q km s '. This value is consistent with the esti- 
mat es of lOu hatha kurta et al. ( 1988 ) based on H i observations 
and lRubinet al.l(ll999l) ^ based on optical observations. 

The main differences between the observed and the mod- 
eled channel maps occur at high velocities, with the model be- 
ing unable to reproduce the emission closest to the center seen 
up to ±150kms Moreover the model cannot reproduce the 
continuous shift of the emission to the northwest and therefore 
the central velocity gradient present in the velocity field. 

In addition, to obtain emission that is spatially elongated 
up to high velocities, we must set a velocity dispersion of at 
least 40kms"' in the ring. This is a high velocity dispersion, a 
fraction of this may come from a shock, but beam averaging of 
the velocities at the ends of the orbits can also be expected to 
contribute significantly. This consideration implies that the real 
orbits must have sharper ends than ellipses. 

3.4.3. Interpretation 

The orbit pattern and the face-on view of the emission strongly 
resemble gas kinematics in a barred potential as observed in 
many galaxies. The large gradient of the orbit position angle 
at 600 pc mimics a tightly wounded spiral forming a ring. The 
change in sign of the gradient of the position angle is a re- 
markable feature that is known to occur at the Inner Lindblad 
Resonance (ILR) when there is only one, or betw een ILRs 
when there are two of them (see e.g. Wadal 1994 ). The fact 
that this feature coincides with a column density peak, which 
is expected at resonances, indicates that this is actually a sin- 
gle ILR. The two symmetric peaks in the face-on integrated 
map and the large velocity dispersion at the resonance would 
then correspond to the shocks frequently observed at the lo- 
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cations of dust lanes in bars, which are due to the crowdin g 
of orbits that favors cloud collisions (see e.g. ICombesI Il996l) . 
Identification of the ILR is further supported by the orienta- 
tion of the bar inferred from the model, which is similar to that 
of the large scale stellar bar. Indeed, th e shocks are expected 
to lead the bar (see e.g. Combes 1996^. which would there- 
fore have a position angle somewhere between and 30deg 
(see Fig. [8] and see also the maximum value for the orbit po- 
sition angle in Table|2l 210deg=30 mod 180 deg). This is con- 
sistent with the large scale stellar bar position angle of 15 deg 
jjogee et al.ll2005b . This agreement implies that the resonance 
observed is unlikely to result from the action of an inner bar, 
but is rather the ILR of the large scale bar. 

Is this interpretation supported by a model based on gas 
orbits in a barred potential? Why is the velocity field with the 
strong velocity gradient twist in the center not reproduced? Is 
this due to the elliptical approximation of the orbits or are pure 
radial motions (as opposed to radial motions intrinsic to non- 
circular but closed orbits) also contributing? To answer these 
questions, we compare the observations to a model based on 
analytical solutions for the gas orbits in a barred potential. 

3.5. Barred potential model 
3.5.1. Description 

This m odel was int rod uced and develop ed by 

Sanders & Huntl evl (Il975h. iLindblad & LindbladI d 1 994 . 
Wadal (ll994 . ISakamotoetal] ( Il999l) and lBakej ( I2OO0I) . It is 
based on analytical solutions for the orbits of gas particles 
obtained by a perturbation calculation of the epicyclic stellar 
motions. The potential is described by an axisymmetric 
component and a weak m-2 perturbation. The perturbation 
has the same radial profile as the axisymmetric component 
and is modulated by a sine in azimuth. The amplitude of 
the perturbation relative to the axisymmetric component is 
defined as the strength of the bar. For the radial profile of the 
potential we use the logarithmic shape. This choice is justified 
by the ratio of the resonance radii. Indeed, the corotation 
radius of the large scale bar, assu med to be 1.2 times t he bar 
extension, is ~5kpc according to Laurikainen & Salol (2002) 
and according to the elliptical orbit model the ILR radius is 
600pc, implying /?c/^ilr~8. As shown in AppendixiBl for 
a logarithmic potential this ratio is 6.3. For co mparison a 
Plummer potential gives a ratio of ~3 ( WadallT994 ). 

As in iBakeJ toOO), two dissipation terms are introduced 
that characterize the damping of radial and azimuthal oscilla- 
tions and are assumed to reproduce the dissipative behavior of 
gas particles. These terms cause the major axis of the orbits to 
lead the bar, the offset being maximal at the ILR. The limita- 
tions of the model are (1) it is singular at corotation; (2) the 
epicyclic motions used to compute the perturbations are based 
on the axisymmetric component only so that the model is not 
fully self-consistent; and (3) computing closed orbits assuming 
dissipation rates is not self-consistent, because dissipative par- 
ticles imply open orbits. As a consequence, the radii of the reso- 
nances obtained with this model (AppendixlBli cannot be taken 
as accurate results. However, this can be seen as a first order ap- 



proximation that works only when the real orbits are near to be- 
ing closed. In this case the resonance radii obtained give a basic 
picture of the scales involved. The orbit patterns produced by 
this model exhibit some of the main features expected in barred 
potentials, e.g., the elongated orbits, the phase shift change at 
an ILR, and the spiral pattern. Also, they match observations 
and r esults of numerical simulations (see Lindblad & Lindbladj 
Il994i) . This analytical approach should be particularly suitable 
for NGC 4569 which, as noted in Section|2] is rather symmetric 
and dominated by an m=2 component in CO and shows a weak 
stellar bai- in the NIR jLaurikainen & Salol |2002|) . The results 
of the elliptical orbit fit allow us to fix some of the parame- 
ters describing the potential. More details about this model are 
given in Appendix lA.4l 

3.5.2. Results 

The best fit parameters for 7 radii, including the parameters 
fixed by the results of the elliptical orbit fit, are listed in Table[3] 
Some of the channel maps of the simulated data are shown in 
the last column of Fig.|6l the integrated map and velocity field 
are shown in Fig.|2]d. The face-on view of the model solution 
is shown in Fig.[8]c and the orbital pattern in Fig.|9l The best fit 
bar position angle of 20 deg is consistent with the large scal e 
stellar position angle of 15 deg derived bv lJogee et al.l (l2005b . 



This result confirms the conclusion drawn from the elliptical 
orbit model, namely that the kinematics inside the inner kilo- 
parsec are driven by the gravitational potential of the large scale 
stellar bafl 

For this model as for the others, a column density jump 
of an order of magnitude is required at a radius of 600 pc and 
a central concentration of similar column density is required. 
A velocity dispersion of 30 km s"' is required at the radius of 
the column density increase. Although this is still a high value, 
it is less than for the elliptical orbit model because the orbits 
have much sharper ends (Fig.|9]l that contribute significantly to 
the velocity dispersion through beam averaging. This suggests 
that the orbit shapes in this model are more realistic than the 
ellipses. 

As can be seen in Fig.|6] this model reproduces the data bet- 
ter than the elliptical orbit model, especially at high velocities. 
The low level emission is not so well reproduced but this is 
expected as this model diverges at corotation, making it unreli- 
able at large radii. The velocity field (Fig.|7]i is also the closest 
to the data when compared with the two other models discussed 
above. 



3.5.3. Interpretation 

The fact that an analytical model that includes certain physical 
assumptions (about the potential shape and gas properties) re- 
produces the data well lends support to its underlying assump- 
tions. In particular, the interpretation of the column density in- 



' Note that although this conclusion drawn from the elliptical orbit 
model was used as an hypothesis to fix some of the parameters of the 
model, these fixed parameters are not related to the bar position angle, 
so the argument is not circular. 
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Table 3. Barred potential model parameters 



galaxy position angle (deg) 20 

characteristic length (pc) 400° 

characteristic velocity (km s"') 220" 
bar position angle (deg) 20 

bar strength 0.015 
pattern speed (km s"' kpc"') 60" 

radial dissipation rate (km s"' pc"') 0.002 
azimuthal dissipation rate (km s"' pc"') 0.005 



charact. radius (pc) 


10 


250 


400 


600 


700 


1500 


3000 


col. dens, (no units) 


1.5 


0.1 


0.5 


1.0 


0.05 


0.02 


0.02 


velocity disp. (km s"') 


35 


10 


30 


30 


10 


10 


10 


scale height (pc) 


10 


10 


10 


10 


10 


10 


10 



value fixed by the elliptical orbit model solution 



crease at 600 pc as an ILR of the large scale bar is well con- 
firmed by this fit. This also supports the assumption that there 
is only one ILR or that the ILR region is small with little space 
for X2 orbits. If there were two ILRs, the gradient of the orbit 
position angle would change in sign at a similar radius (in the 
middle of the 2 ILRs as described in IWadall994 ). but the max- 
imum position angle of the orbits would be much greater, lead- 
ing to a much greater position angle for the oval ring unlikely 
to fit the observations equally well. The good fit also shows 
that even though this model is not fully self-consistent, the or- 
bit pattern and orbit shapes favored by its best-fitting version 
may be reasonably close to reality. We note that the configura- 
tion of the orbit pattern implies that the velocity field is very 
sensitive to the orbit shapes. Indeed, the most elongated orbit 
are almost aligned with the major axis (the bar is almost aligned 
with the major axis), and the strongest velocity components are 
therefore orthogonal to the line of sight. As a consequence the 
velocity field is not only dominated by the elongation of the 
orbits but it is also sensitive to the orbit shapes, i.e. how the ve- 
locity vectors behave at the ends of the orbits. This can explain 
why the diffence between the ellipses and these more realistic 
orbits implies a significant difference in the velocity field. 

4. Mass inflow 

Two mechanisms contribute to angular momentum dissipation 
and mass inflow: gravitational torques, which depend on the 
potential and on the individual gas orbit shapes and orienta- 
tions, and viscosity torques, which depend on the cloud-cloud 
collision rate and on the interstellar turbulence. In the previous 
section the observed molecular gas kinematics were success- 
fully modeled by analytical solutions for gas particle orbits in 
a barred potential. We now address the question of the mass 
inflow based on this model solution. 

As the potential of the model solution is known at each 
point in the galaxy, the gravitational torque can be directly 
computed for each particle following Eq. |C.2| In addition, as 
the orbits are known, the torques can be integrated over the 
orbits to estimate the angular momentum loss over a rotation 
(Eq. lC.4l i as well as the corresponding radial velocity of the 
matter inflow (Eq. lC.6l ). As already outlined, this approach is 



not fully self-consistent, since the orbits should not be closed if 
the angular momentum changes over a rotation, but the angu- 
lar momentum loss thus computed can be considered as a first 
order approximation to the case where orbits are nearly closed, 
i.e., the angular momentum loss due to the gravitational torque 
is a small fraction of the angular momentum. The good fit of the 
model to the observations suggests that, on average, the molec- 
ular gas in the inner kiloparsec lies in this regime. 

The torque is proportional to the bar strength, which is only 
poorly constrained by the 'weak bar' regime requirement (see 
discussion above), to e<0.05. This implies that only upper lim- 
its can be derived for the mass inflow. Figure[TO] shows the net 
torque (as defined Eq. lC.2l ) for e=0.05 on the top left, and the 
corresponding radial velocity on the top right. Both curves look 
similar, because the radial change of the angular momentum as 
determined by the axisymmetric potential is almost constant 
above 500 pc (see Eq. lC.6l ). Due to their circular shape, the or- 
bits with a characteristic radius less than 300 pc undergo no 
net torque and as a consequence the large scale bar would not 
be able to fuel the nuclear starburst on these scales via grav- 
itational torques only. The maximum radial velocity is of or- 
der lOkms"' and occurs for the orbit with a characteristic ra- 
dius of 650 pc, i.e. at the ILR. The value of this upper limit 
for the radial velocities is consistent with our working hypoth- 
esis, namely, that to first order, radial motions can be neglected 
compared to rota tional motions, consistent with the results of 



Wong et alJ(l2004l) . The variations of the radial velocity, shown 



on the bottom left of Fig.[TOl imply that the radial flow is not 
steady over the disk and as a consequence that matter must ac- 
cumulate in some orbits and be depleted in others. The actual 
mass flow also depends on the gas mass available at each ra- 
dius (the radial velocity would directly translate into mass flow 
only if the gas mass did not depend on the radius, i.e., for a 
column density profile following '). Nevertheless, without 
any assumption about the radial density profile, the variation 
of the radial velocity shows that the gas will tend to be de- 
pleted just outside the ILR (the inward velocity increases with 
decreasing radius down to 600 pc) and to accumulate inside the 
ILR (the inward velocity decreases with decreasing radius in- 
side 600 pc). This model is therefore consistent with the for- 
mation of rings at ILR as observed in many spiral galaxies 
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Fig. 10. a) Net gravity torque, b) radial velocity, c) radial velocity gradient and d) mass accumulation rate versus characteristic 
radius of the orbits (see Section|4]for details). All the quantities computed for the maximum bar strength, e=0.05. 



dButa & Combes Il996h and often used as a signature of the 
ILR. 

The mass accumulation rate per radial element at a given 
radius is the gradient of the radial momentum (product of the 
radial linear density with the radial velocity); it is plotted as 
a function of the characteristic radius for e=0.05 in Fig.fTOld. 
To obtain the correct mass units, the column density profile 
of the model was normalized so that the total mass of gas 
equals that measured in Section|2] We see that the maximum 
accumulation rate inside the ILR is of order 0.2MQyr"' pc"'- 
Outside and close to the ILR the depletion rate is of similar am- 
plitude. The accumulation rate integrated inside the ILR (i.e. 
over the radius range 300-600 pc) gives an inflow rate of order 
20Moyr This is significantly higher than other inflow rates 
estimates in the literature which are in the range 1-4 Mp yr ' 
dOuillen et al.lll995l:lLaine et al.lll999l:lRegan et al.lll997|[) . But 



this is an upper limit for this specific model and lower val- 
ues are not excluded. The accumulation of molecular gas in- 
evitably leads to a density increase and star formation. If the 
fraction of infalling gas transformed into stars (star formation 
efficiency) is constant and is of order 10%, then the star forma- 
tion rate (SFR) integrated over the accumulation region would 
be of the order of a few solar masses per year, which is realistic 
when compared to what is observed at ILRs in other galax- 



ies. In NGC4569 however, the Ha morphology and the CO(2- 
l)/(l-0) line ratio suggest that most of the star formation ob- 
served in this galaxy is taking place inside the inner lOOpc (see 
Section|231l. The presence of ~ 5xlO"* O and B stars in the in- 



ner 30 pc (iGabel & Bruhweiler 2002) would imply a constant 



SFR of order 2M0yr (assuming a Salpeter IMF and a life 
time of 5 Myr for O and B stars). This implies either that an- 
other fueling mechanism is taking over inside the inner 300 pc, 
or that the nuclear starburst is ending and star formation at the 
ILR will dominate in the future. 

Viscosity torques are intrinsic to the model and are spec- 
ified through the radial and azimuth al dissipative parameters. 
However, as outlined in lBaker ( 2000 V these parameters lead to 
closed orbit solutions, so they cannot be directly used to com- 
pute the dissipation of angular momentum and the mass inflow 
rate. Estimating the intensity of the viscosity torques based on 
the density and velocity dispersion at a given point in the disk 
would require adopting prescriptions for the visosity that are 
generally acknowledged as very uncertain. Thus, our model 
does not allow us to directly constrain the rate of mass inflow 
due to viscosity torques. It could well be that these torques take 
over the fueling of the nucleus inside the 300 pc where the grav- 
ity torques are not efficient anymore. Alternatively, an inner 
stellar bar could drive the gas further inward as in NGC 6946 
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dSchinnerer et al. I l2006h : however, the current observations do 
not show any evidence of such a small inner bar. 

5. Discussion 

Models of gas dynamics in galaxies come in two main cat- 
egories; ballistic models in which th e gas is represented b 



egories; ballistic moaeis in wnicn tn e gas is representea Dy 
colhding clouds on ballistic orbits ( Combes & GerinI 
Bvrd et al .11994 . and hydrodynamic models in which the gas is 
assumed to be a c ontinuous dissipative medium (Athanassoula 
1992t IPiner et alJ [19951: IWada & Habd Il995l: iFriedli & Ben j 



1993 '). Comparisons of the two approaches are discu s sed in 
yndblad & Lindblad (1994), Wada (1991, ICombesI (|l996|) 
and [Regan et al.l ( Il999h . The analytical model of orbits in a 
barred potential presented in the previous section effectively 
traces the evolution of gas particles in the gravitational poten- 
tial as they are subjected to collisions. Although the equations 
are continuous, they are not based on hydrodynamic assump- 
tions, and the gas is assumed to be an ensemble of independent 
particles rather than a continuous medium, meaning the model 
can be classified as ballistic. The good fit of this model to our 
CO observations presented here seems to validate such a bal- 
listic approach. The results of the two approaches may actually 
not have so many differences t hat the observati ons can distin- 
guish. As already suggested bv lLindblad & Lin dblad ( 1994), it 
may be that the orbits out of this analytical model describe a gas 
flow similar to the one obtained with a hydrodynamic model. 
This hypothesis is supported by comparing the orbit pattern 
as shown in Fig.|9]b to the gas flow resu lting from a hydro- 
dynamic simulation as shown in Fig. 9 of iRegan et al. ( 1999 ). 
The shapes of the orbits resemble the gas streamlines in their 
model and the dust lane/shock region of their model is well 
reproduced by the orbit crowding in our model. The main dif- 
ference is that the largest circular orbit in our model is located 
well inside the ILR, whereas it defines the ILR in their model. 
As a consequence, the ring formed at the ILR by gas accumu- 
lation is elongated in our model (Fig. [8lc) rather than c i rcular . 
But such a shape has also been seen by lNakanishi et al.l (120051) 
based on Smoothed Particle Hydrodynamics (SPH) simulations 
of the same galaxy. 

The differ ence between ballis tic and hydrodynamic mod- 
els outlined in Regan et al. ( 19991) corresponds to the sharp in- 
version in the radial component of the velocity vector in the 
shock/dust lane region. This feature is argued to be reproduced 
by the hydrodynamic models but not by the orbits of ballis- 
tic models. However, this argument is based on the assumption 
that the orbits are elliptical (see their Fig. 10). As emphasized in 
Section [33] the analytical orbits of gas clouds in a barred poten- 
tial used here have much sharper ends than ellipses, blurring the 
nominal distinction between hydrodynamic and ballistic mod- 
els. 

It is also interesting to note that the conclusions of both 
approaches in terms of mass inflow are not so different. The 
inflow stops at the largest circula r orbits (the "nuclear ring") in 
the model of iRegan et al. d 19991) an d at the largest x? o rbit or 
at the largest nonlooping Xi orbit in Regan et al.l (120041) . In all 
cases, the bars seem to be incapable of driving molecular gas 
inside a radius of 100 pc. Similarly, in our model, the gravita- 



tional torques stops where the orbits become circular at a radius 
of 300pc. 

The advantage of the method followed here is that we are 
able to first identify the ILR based on a kinematic argument, 
namely the inversion of the radial variation of the ellipse orien- 
tations in the second model (Section lSTil i. This allows us to fix 
the potential parameters of the third model, i.e. the orbit pat- 
tern, and to fit the radial profile of the column density and ve- 
locity dispersion independently. Although limited in its domain 
of application (weak bars and inside corotation), this method 
has the advantage of being insensitive to initial conditions or 
numerical effects that may affect evolutionary models (ballis- 
tic or hydrodynamic). In addition, it is obviously very unde- 
manding in terms of computation resources and well adapted 
to iterative fitting to observations. 

6. Conclusion 

The main results of this article can be summarized as follows: 

1. In NGC 4569, 70% of the [1.1 +0.2] x IO'^Mq of moleculai- 
gas detected in the inner 20" is concentrated within a radius 
of 800 pc and distributed along the large scale bar, with two 
peaks at ~500pc from the center and one peak close to the 
center. There is a hole in the CO distribution at the location 
of the nuclear starburst where evidence of outflow is seen 
in Ha emission. 

2. Fitting a model based on elliptical orbits reveals a change in 
sign of the radial variation of the orbit orientations at a ra- 
dius of ~600pc, where molecular gas is also accumulated. 
This coincidence, together with the value of the orbit posi- 
tion angle maximum, leads us to identify this feature with 
confidence as a single ILR of the large scale bar. 

3. The kinematics are well reproduced by a model based on 
analytical solutions for the gas particle orbits in a weakly 
barred potential with a single ILR at 600 pc. 

4. The gravitational torques implied by this model are able to 
efficiently funnel the gas inside the ILR down to a radius 
of 300 pc but another mechanism must take over to fuel the 
nuclear starburst, i.e., there is no evidence for ongoing fu- 
eling of the starburst by the large scale bar If such a mech- 
anism does not exist, the nuclear starburst will stop and star 
formation at the ILR will dominate with a star formation 
rate of order one solar mass per year (assuming a 10% star 
formation efficiency). 



This last point agrees with the result of lGarcia-Burillo et al 



(|2005|) , namely that the gravitational torque caused by kilopar- 
sec scale perturbations of the potential seems not to be the only 
mechanism responsible for fueling galactic nuclei. 

Appendix A: Detailed description of the modeling 

A.1. All models 

The first two models are entirely specified by a radially vary- 
ing set of pr operties similar to th at of the widely used 'tilted 
ring model' (^Rogst ad et al. 1974 ). The values of the parame- 
ters are specified at some galactocentric radii and linearly in- 
terpolated between those radii. The third model has additional 



F. Boone et al.: Molecular gas in NUclei of GAlaxies (NUGA) 



15 



global parameters. Thus, the total number of free parameters 
is adjustable through the number of radii, A^iadii, and is given 
by A^radii X A^properties + A^giobah where A^propeities IS the number of 
properties specified at each radius including the value of the 
radius. The number of radii required to reproduce the data sat- 
isfactorily was found to be greater than 7. In this case, for all 
models the number of free parameters is too high compared 
to the dynamic range of the data (signal-to-noise ratio, spatial 
dynamic range and spectral dynamic range) to permit a rigor- 
ous minimization of an objective function like over the full 
parameter space. Since the goal of this modeling is to test the 
compatibility of the different hypotheses with the main features 
of the data, rather than to estimate to a high accuracy the nu- 
merous free parameters, the comparison of the simulated data 
to the real data has been done by ey^l 

In all the models, the inclination of the galaxy is 
fixed to 70deg JtuIIv 1981, its position angle to 30deg 
dGuhathakurta et al.nl 988i, except in the third model where 
the PA is also a free p arameter) and its distance to 17Mpc 



dBottinelli et al.l 119831: iTuUvi 119881 only used to convert an 



gular scales into physical scales). In addition, in all models we 
make the eff'ective assumption that mass column density is pro- 
portional to the CO intensity with the same constant of propor- 
tionality everywhere. 

A.2. Axisymmetric model 

In the first model, the properties specified are the radii, the ro- 
tational velocity, the radial velocity (positive in outward direc- 
tion), the column density through the disk, the disk scale height 
(FWHM of the Gaussian profile) and the local velocity disper- 
sion (square root of the Gaussian distribution variance). The 
number of free parameters is therefore 6 x A^radii- The vertical 
distribution of matter and the velocity distribution at a given 
position are assumed to be Gaussian. 

A.3. Elliptical orbit model 

To mimic the formation/destruction of the clouds we apply an 
exponential taper to the column density along the orbits. The 
taper equals 1 at the aphelion and decreases like exp(-20/7r), 
where is the azimuth along the ellipse. The sign of 6 is de- 
termined by the sense of rotation of the cloud in the reference 
frame of the density wave for which a pattern speed needs to be 
specified. Here we assume that all the gas lies within the coro- 
tation radius, and that the density decreases smoothly in the 
sense of rotation of the whole galaxy, i.e. counter-clockwise. 

The properties of the model are the characteristic radius, the 
characteristic rotation velocity (the angular momentum), the el- 
lipticity, the position angle of the ellipse, the column density, 
the velocity dipersion and the scale height. The pattern speed of 



^ This comparison uses the DALIA interface to numerical simula- 
tion codes which is being developed jointly by the Observat oire de 
Paris and the Max Planck Institut Fiir Radioastronomie (Boon e^t al.| 
l2006h . This graphical user interface gives access to the model param- 
eters and allows the user to directly compare the outputs of the code 
with real observations. 



the entire model was not well constrained by the observations 
and since the larger the value the worse the fit, it was decided 
to fix this parameter at zero. The number of free parameters is 
therefore 7 xA^radii- As in the previous model, the vertical distri- 
bution of matter and the velocity distribution at a given position 
are assumed to be Gaussian. 



A.4. Barred potential model 



We use 2 parameters describing the logarithmic potential (the 
characteristic scale and velocity, see Eq. lB.lb . 3 parameters for 
the bar peiturbation (the position angle relative to the galaxy 
major axis, the strength, and the pattern speed) and 2 parame- 
ters for the gas dissipative properties (the radial and azimuthal 
dissipation rates). We found it difficult to obtain a good solu- 
tion with a galaxy position angle of 30 deg as in the previous 
models, so we left this as an additional free parameter here. In 
addition, the gaseous disk is described by 3 properties (the col- 
umn density, the velocity dispersion, and the scale height) at 
different radii as in the previous models (the radius here is the 
characteristic radius of the orbit, which coiTesponds to the ra- 
dius of the circular orbit without bar perturbation). The number 
of free parameters is therefore 8 -n 4 x M-adii- 



The results of the elliptical orbit fit allow us to fix some of 
the parameters describing the potential. First, the characteristic 
velocity of the axisymmetric component of the potential corre- 
sponds to the rotation velocity at infinity. Based on previous re- 
sults, we take v^- 220 km s"' , also consistent with the value of 



Guhathakurta et al. (11988 ') based on Hi observations. Next, the 
fact that the ILR was identified at 600 pc from the center in the 
elliptical orbit model allows us to fix the characteristic scale of 
the potential (using Eq. |B.5t to = 400 pc. Finally, assuming 
there is only one ILR (as the elliptical orbit model suggests), 
the pattern speed is set by Eq. |B.6| to Qp=60kms ' kpc The 
remaining global parameters (as opposed to those specified at 
each radius) are the bar strength, its position angle and the 2 
dissipation rates. There is a balance between the bar strength 
and the dissipation rates, in the sense that the stronger the 
bar or the smaller the dissipation rates the closer the orbits 
are to pure stellar orbits. Stellar orbit s are organized in fam- 
ilies, the so-called xi and X2 orbits , dContopoulos & Grosbol 
1989t ISellwood & Wilkinsonlll993b . which are orthogonal to 
each other, intersect, and can exhibit large loops. As long-lived 
gas orbits cannot intersect or have loops, a realistic solution for 
the orbits must keep the bar strength and dissipation rates bal- 
anced. The bar strength is required by self-consistency to be 
in the weak perturbation regime (less than 5%) but below this 
threshold there is a degeneracy with the dissipation rates. The 
ratio of the two dissipation rates has little effect on the orbits, 
i.e. there is a degeneracy between them. So these three param- 
eters (bar strength and the two dissipation rates) are not well 
constrained and only their combination matters. 
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Appendix B: Resonance radii and pattern speed in 
a logarithmic potential 



A logarithmic potential is defined as: 
1 

Oo(r)= -y^ln(l + -) 



where Vp and rp are the characteristic velocity and length of the 
potential. The rotation frequency is given by: 



rp + 



, 1 d ((Do) 
r or 

and the epicyclic frequency by: 



K{r) 



dr 



(B.3) 



The inner Lindblad resonance (ILR) occurs when f2 - Kj2 
equals the pattern speed of the bar If we assume only one ILR, 
this happens when Q - /c/2 is at its maximum. The radius of the 
ILR, /?iLR, is therefore obtained by solving 



d(Q.-KlT) 



0. 



(B.4) 



It scales with the scale radius of the potential, rp, and is given 
by 



ri=/?iLR/rp=5'/4 



(B.5) 



where r\ is the ILR radius in units of rp. 

The pattern speed of the bar, Qp, is given by Qp - 0(/?ilr)- 



1 ^^i<iT'\ 



— ^ 0.11 — 



(B.6) 



The corotation resonance occurs at the radius where the ro- 
tational frequency equals the pattern speed of the bar. So, the 
corotation radius, R^, is related to R\yr. through 



^(R,) = 0(/?ilr) - ^(i?ILR)/2, 

which leads to: 



rl = 6rf +19+22rr2+8^-4^(4^g^-2^4^-4^ ^2rf + 6rf + 4(B.8) 

where - Rdr^ is the corotation radius scaled to rp. 
Replacing ri by Eq. lB.SI gives: 



re = a/20 VS +43. (B.9) 
The ratio of the resonance radii is given by Eq. lB.5l and |B.9l 



RJR 



ILR 



a/20V5 



+ 43 



51/4 



^ 6.26 



(B.IO) 



Appendix C: Mass inflow in a logarithmic potential 
with a bar perturbation 

The non-axisymetric component of the potential (the bar com- 
ponent) is given by 



(B.l) (Db(r, 0) = eOo cos(2<;*). 



(C.l) 



where is given by Eq. lB.ll and e is the strength of the bar. 
The resulting gravitational torque exerted on a particle at (r, (p) 
is 



(B.2) T{r, (p) = rF^ = 



50 



--ev^ln 



1 + 



sin(2(;A) 



(C.2) 



We define the net torque, 7~, as the average torque over an orbit 
llT{t)dt 



T 



(C.3) 



where t is the orbital period. The angular momentum variation 
over this orbit and per unit mass is then given by 



AM = Tt. 



(C.4) 



A particle that changes angular momentum will settle on a new 
orbit which has an intrinsic angular momentum (as determined 
by the axisymmetrical component of the potential) correspond- 
ing to this new value. Hence, the radial velocity of the matter 
flow is given by the ratio of this angular momentum variation to 
the one fixed by the axisymmetric component of the potential 
and divided by the rotation period. 



Vi-ad = 



AM 1 



T 



dMIdn dMIdr' 
where dAi/dr - d(r^Q.)/dr. Using Eq, 
t 

Viad 



this leads to 



rn(2 - r2/(r2 + r^)) 



(C.5) 



(C.6) 
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